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1. Introduction 



1.1. Orientation 

Fusion rules yield the number of independent couplings between three given primary- 
fields in conformal field theories. We are interested in fusion rules in unitary conformal field 
theories that have a Lie group symmetry, that is, those whose generating spectrum algebra 
is an afiine Lie algebra at integer level. These are the Wess-Zumino-Witten models [1,2]. 
Primary fields in these cases are in 1-1 correspondence with the integrable representations 

(k) 

of the appropriate affine Lie algebra at level k. Denote this set by and a primary field 
by the corresponding affine weight A. Fusion coefficients ^ are defined by the product 

AxA= E M^^P (1.1) 



(For a review of conformal field theory and in particular fusion rules, see [3]; to a large 
extend we follow the notation of this reference.) 

In the infinite-level limit and for fields with finite conformal dimensions, the purely 
affine condition on weight integrability is relaxed and the primary fields are solely charac- 
terised by their finite part, required to be an integrable weight of the corresponding finite 
Lie algebra. Recall that a finite weight A is characterised by its expansion coefficients in 
terms of the fundamental weights coi 

r 

A = ^Aia;i - (Ai,...,A^) (1.2) 

i=l 

where r is the rank of the algebra. The numbers A^'s are the Dynkin labels. The set of 
weights with non-negative Dynkin labels (the integrable weights) is denoted by P+. 

In the infinite-level limit the fusion coefficients reduce to tensor-product coefficients: 

limAAf)^=jVA/. (1.3) 

fc— >oo '^A' 



where N'x^'^ is defined by 



A®/i= J] AAa/jv (1.4) 
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By abuse of notation, we use the same symbol for the highest weight and the highest- weight 
representation. Notice that 

Mx^''=^fx^.* (1.5) 

where z/* denotes the highest weight of the representation conjugate to that of u. Equiv- 
alently, Mxjxu* gives the multiplicity of the scalar representation in the triple product 

A tensor-product generating function codes the information for all the tensor products 
of a given algebra in a single function defined by 

G{L,M,N)= ^f\^.''L^M^'N'' (1.6) 

where = L^^ ■ ■ ■ L^"^ and similarly for and N'^. G can generally be expressed as a 
simple closed function of its variables. For instance, for stt(2), it reads 

^' ^) = (1-LM)(1-LN)(1-MN) f^-^) 

An example of basic global information that can be deduced from a generating func- 
tion is the integrality as well as the positivity of the tensor-product coeflBcients. More 
importantly, from our point of view, is that in the context of fusion rules, the construction 
of the simplest generating functions led to the discovery of the notion of threshold levels 
[4]. Moreover, as shown in the sequel paper, setting up a fusion generating function is a 
way to obtain explicit expressions for these threshold levels. Our new approach to fusion- 
rule generating functions, which originates from the generalisation of techniques developed 
in the present paper on tensor products, leads to a further new concept, that of a fusion 
basis. 

1.2. Overview of the paper 

The present article is organised as follows. We start by explaining in detail the con- 
struction of tensor-product generating functions for finite Lie algebras. The first construc- 
tion which is presented is the character method developed by Sharp and his collaborators 
(section 2). 

Although it is conceptually very simple, the character method is limited by its inherent 
computational difficulties: the disproportion between the simplicity of the resulting form 
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of the generating function and the intermediate calculations is enormous. This motivates 
our alternative approach to the construction of tensor-product generating function. It is 
based on the reformulation of the problem of calculating tensor products in terms of the 
solution of a set of linear and homogeneous Diophantine equations (cf. section 3). The 
elementary solutions of these Diophantine equations represent "elementary couplings" . For 
sp(4), the use of the Berenstein-Zelevinsky inequalities to obtain the elementary couplings 
and their relations (cf. the analysis of section 6) is new. 

The key difficulty is finding the numerous relations that exist in general between the 
elementary solutions. From the Diophantine-equation point of view, the decomposition 
of a solution may not be unique because different sums of elementary solutions could 
yield the same result. To solve this problem we first "exponentiate" it: given a solution 
a. = (cKi, . . . , CKfe) to our system of linear Diophantine equations, we introduce formal 
variables Xi, . . . and consider the monomial Xf^ . . ■X'^^ . The linear span, i?, of all 
such monomials is a "model" for the generating function for the solutions to the original 
set of linear Diophantine equations (see section 5), since the Poincare series of R is the 
required generating function. This series can be calculated using Grobner basis methods. 

For su{N) there is a remarkable graphical construction for computing tensor product 
multiplicities, the famous Berenstein-Zelevinsky triangles. These are introduced in section 
6. We also discuss the analogous construction for sp(4), whose diagrammatic representation 
is new. But the main interest of these re-formulations is that it yields a simple and 
systematic way of obtaining the elementary couplings from the construction of a vector 
basis. Thus we get a new way of constructing the corresponding generating functions. 

2. Generating-function for tensor products: the character method 

2.1. The character method for the construction of the tensor-product generating function: 
the su{2) case 

The method developed by Sharp and collaborators for constructing generating func- 
tions for tensor products is based on manipulations of the character generating functions 
[5]. Although simple in principle, these manipulations become rather cumbersome as the 
rank of the algebra is increased. To illustrate the method, we will work in complete detail 
the simplest example, the su{2) case. 
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The first step is the derivation of the character generating function. The Weyl char- 
acter formula for a general algebra of rank r and a highest-weight representation A is 

XA = ^ (2.1) 

where p is the finite Weyl vector, p = Yli=i '^ii where the characteristic function ^ is 
defined as 

Ca+p = (2.2) 

where e{w) is the signature of the Weyl reflection w and W is the Weyl group. 
For su(2), W contains two elements: With 

X = e'^i (2.3) 

the su{2) characteristic function ^ for the representation of highest weight rmvi = (m) is 

xm+l _ (2.4) 

The character reads then 

Xm = = , 2 = ^"^ + + • • • + ^-"^ (2-5) 

X — X ^ 1 — X 

The character generating function xl is obtained by multiplying the above expression by 
L"^ where L is a dummy variable, and summing over all positive values of m: 

oo ^ oo 

Xl{x) = V L^Xm = =1 V L^{x^+' - x-^-') 

X — X ^ ^—^ 

(2.6) 

1/1 x-^ 



\-x-'^\\-Lx l-Lx-^J {I - Lx){l - Lx-^) 

We should point out here that in all generating finictions in this paper, expressions of 
the form 1/(1 — a) should be formally expanded in positive powers of a. So for example, 

1/(1 — Lx~^) = 1-1- Lx~^ + L'^x~'^ + By construction, the character of the highest 

weight (m) can be recovered from the power expansion of xl as the coefficient of the term 
L"^. The characteristic generating function is defined by 

XL{x) = ^f (2.7) 



and it reads 



^^^^^ {\-Lx){\-Lx-^) 1-Lx 1-Lx-^ ^^"^^ 
the last form being the one that results directly from (2.5). 

The tensor product of two highest-weight representations can be obtained from the 
product of the corresponding characters: 

XmXn = ^^fmJ Xt (2-9) 

I 

This information can be extracted from the product of the corresponding generating func- 
tions. We are thus led to consider the product Xl{x)xm{x). To simplify the analysis of the 
resulting expression, notice that the information concerning the representations occurring 
in the tensor product is coded in the leading term of the character, i.e., the term x"^^^ . 
To insure that every positive power of x singles out a highest-weight representation, we 
can multiply both sides by ^o- To read off these terms, we can focus on the terms with 
strictly positive powers of x in the product Xl(^)Xm(3^)^o(^)- If we require the Dynkin 
label of the representations (and not their shifted value), it is more convenient to divide 
by X before doing the projection, now restricted to the non-negative powers of x. The 
truncation of an expression by its negative powers of x will be denoted by the MacMahon 
symbol [6] f2, defined by 

oo 

f2 J] CnX^ = J] CnX" (2.10) 
— — oo n>0 

When there is no ambiguity concerning the variable in terms of which the projection is 
defined, it is omitted from the Q, symbol. 

We are thus interested in the projection of the following expression 

Xl{x)xm{x)^q{x)x~'^ = Xl{x)^m{x)x~'^ 



1 x-^ \ (2.11) 



(1- Lx)(l- Lx-^) \1- Mx l-Mx-\ 

For these manipulations, we use systematically the following simple identities: 

1 1 f 1 Bx- 

+ 



(l-Ax)(l-Sx-i) {1-AB)\1-Ax 1-Bx-^ 



1 f Ax 1 , /„.„x 

+ TT-ZT (2-12) 



(l-AB) \1-Ax 1-Bx 

1/1 1 

+ z — T - 1 



(l-AB) Vl - via; 1 - Bx''^ 
5 



There are two terms to analyse. The first is 

1 1/1 Lx-^ , 

+ (2-13) 



(1 - Lx){l - Lx-i)(l - Mx) (1 - Lx){l - LM) \1-Mx 1 - Lx'^ 
The first part is not aff'ected by the projection and the second can be written as 
Lx~^ Lx~^ ( Lx 



{1 - Lx){l - LM){1 - Lx-^) {l-LM){l-L^) \1-Lx^ l-Lx-^J ^'^'^^^ 

The second term of this expression contains only negative powers of x and can thus be 
ignored and the first part is unaffected by the projection. We have thus, for the first term 
of (2.11) 

1 1 / I \ 

> (I - Lx){l - Lx-''){1 - Mx) ~ {l-Lx){l-LM) \l-Mx^Y^) 

The projection of the second term of (2.11) is: 

x-^ 

n 



> (1 - Lx){l - Lx-^){1 - Ma;-i) 

x-"^ f 1 Mx-^ 

nv-. 7 TV7-. TTTT ^ + 



> (1 - Lx-'^){1 - LM) \1 - Lx 1-Mx-'^ 

= n 



> (1 - Lx-^){1 - LM){1 - Lx) 



X ^ ( Lx 

= g (l-LM)(l-L^) (l^I^+l^i^.l (^-^^^ 

Lx-^ 

= Q 



> (1-LM)(1-L2)(l-Lx) 
Lx-^ ( 1 



(1 - LM)(1 - L2) \\-Lx~^ 

_ L2 

~ {1 - LM){\ - L'^)(\ - Lx) 

Subtracting (2.16) from (2.15), we find that 

^XL{x)iM{x)x-^ = ^^_LM){1-Lx){l-Mx) ^'-'^^ 
Replacing x by AT, we thus get 

^""^'(^■^■^)= (l-LM)(l-LiV)(l-M^) (^-^^J 
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2.2. The abstract setting: Poincare series, elementary couplings and relations; defining a 
model 

As we shall see it is frequently useful have a model., R, for a generating function 
G{Xi, . . . , Xk) such as (2.18). By this we mean a commutative Q-algebra with an identity, 
graded by N'^, ( N = {0, 1, 2, 3, . . .}) 

R = ©aeN^-^a ) RaRp Q Ra+P (2.19) 

and such that its Poincare series (also frequently called Hilbert series) 

F{R) = dimQ(i?a)X" 

satisfies 

F(i?) = G(Xi,...,Xfc). (2.20) 

For example, for (2.18), with Xi = L, X2 = M, X3 = AT, we can take R = Q[Ei, E2, -E3], 
which is the polynomial ring generated by the formal variables Ei,E2, -E3 (in fact all our 
examples R is either a subring or quotient of a polynomial ring) with the grading of Ei , E2 
and E^ being (1, 1, 0), (1, 0, 1) and (0, 1, 1). The homogeneous subspaces are spanned by 
E^E^E^, a,b,ceN with grade {a + b,a + c,b + c) and so 

F{R)= X^+''X^+''X^+'' = G'''^^\X^,X2,Xs) (2.21) 

(a,6,c)eN3 

as required. 

If R is generated by elements Ei, . . .,Eg and is a model for a generating function 
G for tensor products (or fusion products) then we call ^^i, . . . , E'g a set of "elementary 
couplings" for G. 

It should perhaps be stressed that a priori the variables Xi, . . . , Xk and . . . , i?s 
are unrelated. We shall refer to the -E's as model variables and the X's as grading variables. 
If the grading vector of is a* , z = 1, . . . , s then there is an associated monomial in the 
grading variables: X" , for which wc will use the notation g{Ej). For example in the above 
example we have g{Ei) = = LM. However, to avoid tedious repetition when 

writing down generating functions we shall often write, for example, 1/(1 — £■!) rather than 
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1/(1 — g{Ei)). In all such cases where model variables appear in a generating function 
they should be replaced by the corresponding monomial in the grading variables. 

In the case of tensor products we use the notation "i? : g{E) : product" to denote 
a set of elementary couplings with their "exponentiated" grading and the corresponding 
term in the tensor product. So in the example above we would write: 

E^:LM: (1) (g) (1) D (0), 

E2:LN: (1) ® (0) D (1), (2.22) 
Es : MN : (0) ® (1) D (1) 

Having made the distinction between grading and model variables, it should be noted 
that there are cases where we can identify the model as a ring generated by monomials in 
the grading variables. So in the above example we could define Ei = LM, E2 = LN and 
i?3 = MN and take the model for our generating function to be the subring of Q[L, M, A^] 
generated by E2 and £"3. However, it is not always desirable, or even possible, to make 
this identification. 

We close this section with two examples of how models for the su{2) character gener- 
ating function can be constructed. 

The first method, which has been exploited by Sharp et al (see [5]) to construct 
character generating functions, amounts to finding an algebra R which is a module for the 
Lie algebra su{2) and such that, as an ^^(2) module, R is isomorphic to ©i>iVi where Vi 
is the irreducible su{2) module of dimension i. 

In this case we can take R = Q\p, q] with the generators of su{2) being given by 
differential operators: 

d d d d , , 

h=p-^ -q-^, x_ = q—, x+=p— (2.23) 
op oq op oq 

The su{2) highest- weight vectors are p*, i > and a basis of the irreducible submodule 
of dimension i is just given by the monomials of degree i m. p and q. We can give R an 
grading by taking the degree of p to be (1, 1, 0) and of q to be (1, 0, 1). Here the first 
grading index specifies the representation while the other two refer to a particular weight. 
As i? = Q[p, q] the Poincare function for R is, 

1 

{l-p){l-q) 



with the understanding, as explained above, that p and q should be replaced by the corre- 
sponding expression in terms of the grading variables. Let us denote these grading variables 
here by L (which exponentiates the representation index) and x, y (exponentially related 
to the weights). The Poincare function reads then 

(2.24) 



{1-Lx){l-Ly) 

Another way of constructing a model for the weight generating function, which makes 
more natural the grading, is to observe that the complete set SU{2) weight vectors 
of finite dimensional irreducible su{2) modules are in 1-1 correspondence with one-rowed 
Young tableaux. If the Young tableau has c boxes filled with a I's and b 2's then there is 
a constraint 

a + b-c^O, a,b,c>0 (2.25) 

and so the solutions to this linear Diophantine equation are in 1-1 correspondence with 
the complete set of SU (2) weight vectors. Thus to find a model for the weight generating 
function it is sufficient to find a model for the solutions to (2.25). It is not difficult to see 
that every solution to this equation is a linear combination (with non-negative coefficients) 
of the two fundamental solutions: (a, 6, c) = (1, 0, 1) and (a, b, c) = (0, 1, 1). Let R be the 
subring of Q,[A, B, C] generated by the monomials Ei — AC, E2 = BC. Considering the 
exponents of the monomials Ei and E21 we see that the monomials in R correspond to the 
solutions of (2.25) and hence taking the natural grading on R ensures that the Poincare 
series of R is the generating function for the solutions to (2.25) and hence is the required 
generating function. In this example there are no relations between Ei and E2 and so R 
is isomorphic to the polynomial ring in two variables (as expected) and so the Poincare 
function is once again (with A ^ x, B ^ y,C ^ L) given by (2.24). 



2.3. Multiple su(2) tensor products 

In order to illustrate the occurrence of linear relations between elementary couplings, 
consider the problem of finding the multiplicity of a given representation C, in the triple 
product \® jJL® V. In terms of character generating functions, this amounts to consid- 
ering the product Xl{x)xm{x)xn{x) D Xp{x)-, or equivalently, xl{x)xm{x)^n{x)x-'^ D 
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^p{x)x ^. The left side is then projected onto positive powers of x. We are thus led to 
consider 

1 x-^ \ 

^ {1 - Lx){l - Lx-'^){1 - Mx){l - Mx-^) \1 - Nx 1 - iVa;" V ^' ^ 

The projection of each term is worked out as previously and the resulting expression is 
found to be, with x replaced by P: 

G(L, M, N, P) = 1 - LMNP ^ 

^ ' ' ' ^ {1- LP){1- MP){1- NP){1- LM){1- LN){1- MN) ^ ' 

This is the sought for generating function. Here we would like to have a model with 6 
elementary couplings corresponding to the terms in the denominator of the generating 
function: 

Ei-.LM: (1) ® (1) ® (0) D (0) 
E2:LN: (1) ® (0) ® (1) D (0) 
Es-.LP: (1) ® (0) (8) (0) D (1) 

(2.28) 

: MN : (0) (g) (1) (g) (1) D (0) 

E^-.MP: (0) ® (1) ® (0) D (1) 

Ee:NP: (0) ® (0) ® (1) D (1) 
and there must be a linear relation (in this context, such a relation is often called a syzygy 
in the physics literature - see in particular [5] and related works) between the following 
products (signalled by a term in the numerator) which has grading LMNP: 

EiEq, E2E5, E^E^ (2.29) 



It is not difficult to see that a model is given by Q[ei, 62, 63, 64, 65, cq]/! where Ei = 
Cj + /, z = 1, . . . , 6 and / = (aeiee + ^6265 + 06364) is the ideal generated by the polynomial 
a6i66 + 66265 + C6364 for any choice of a, 6, c e Q not all zero. 

The elements of R have the form m + 1 with m G Q[6i, ...,60]. However there is no 
canonical way of choosing the representatives m. Take for example the case a = b = c = 1. 
(Usually we will construct a model for our generating function as explained above and this 
construction will fix the values of a, b and c). In R we have E-^Eq = —{E2E5 + E^E^) and 
so we can take as a basis for R the set of (equivalences classes of ) monomials which do 
not contain the product EiEq. In this case we say that we have chosen to make EiEq a 
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'forbidden product'. Similarly we can forbid the products E2E5 or As we shall see 

later, the choice of forbidden products corresponds to a choice of term ordering. 

Before leaving this example, we would like to rework it from a different point of view, 
as an illustration of the 'composition' technique of generating functions. Let G{L, M, R) 
describe the tensor product corresponding to XlXm ^ Xr similarly let G{Q,N, P) 
correspond to XqXn ^ Xp- We are interested the product Xl{x)xm{x)xn{x) D Xp(^); 
but treated from the product of the two generating functions G. We thus want to enforce 
the constraint R = Q in the product G{L, M, R)G{Q, N, P). The idea - which is used in 
the references in [5] mainly in relation with the construction of generating functions for 
branching functions - is to multiply this product by {1 — Q~^R~^)~^ and, in the expansion 
in powers of R and Q, keep only terms of order zero in both variables: with an obvious 
notation we have 

h?lG{L,M,R)G{Q,N,P) ^ 



= |g 5]An(i^,M)i?'^ J]i?^(iv,p)g-5^i?-V' (2.30) 

n m H 

^Y.Ap{L,M)Bp{L,M) 
p 

which is manifestly equivalent to considering 

hG{L,M,x)G{x-^,N,P) (2.31) 

With the explicit expressions for the generating functions, we have thus 
X I 1 

{1 - Lx){l - Mx)(l - LM) (1 - Px-^){1 - Nx-^){1 - NP) ^^"^^^ 
A brief and by now standard analysis yields directly the generating function (2.27). 

24. The sp{A) case 

As a final example, consider the sp{A) case. With the = z = 1, 2, the charac- 
teristic function is found to be 

t, , _ ^"i+l^T^+l m+n+2 „n+m+5 -n-1 , m+2n+3 -m-n-2 

S(m,n) — "^2 -^1 •*'2 "^1 •*'2 ~r J-i X2 

(2.33) 

I ™-m-2n-3 n+m+2 _ m+1 -m-n-2 _ -m-2n-3 n , -m-1 -n-1 
~r uy -j^ "^2 J-" ^ "^2 t*-" ]^ "^2 ~'~ "2 
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and the characteristic generating function is 

^ 1 

^' ^ {1 - LiXi){l - LiXiX2^){l - L2X2^){1 - L2X^^X2) 

I + L2 , {l + L2)LiXi 



X 



+ n T^TTn TTIT^ (2-34) 



(1 - L2xlx2'){l - L2X2') (1 - Lixi)(l - L2XIX2') 

LiXi^X2 



(1 - Lia;i)(l - Li^i ^a;2). 
From this we construct the character generating function and then we can proceed to the 
tensor-product generating function. This is again extremely cumbersome. The result is [7] 

G'«f(4)(Li,L2,Mi,M2,iVi,iV2) 

= [(1 - MiNi){l - LiiVi)(l - LiMi)(l - M2N2){1 - L2N2){1 - L2M2)]-' 



1 L2M2N( 



2 



(1 - L2AfiiVi)(l - L2MIN2) (1 - L2MiiVi)(l - L2M2NI) (2.35) 
LIMIN1N2 L1M2N1 



+ 



{1 - LiM2Ni){l - LIM2N2) {I - LiM2Ni){l - L2M2NI) 
LIM2N2 L1M1N2 



+ 



' {1 - LiMiN2){l - LIM2N2) {1 - L^M^N2){1 - L2MIN2) , 
From this expression, we read off the following list of elementary couplings (recall that the 
first variable is a model variable and then we write the corresponding monomial in the 
grading variables): 

Ai-.MiNi, A2:LiNi, A3 : LiMi 

Bi : M27V2, B2 : L2iV2, S3 : L2M2 

(2.36) 

Ci : L2M1N1, C2 : L1M2N1, C3 : L1M1N2 
Di : LIM2N2, D2 : L2MIN2, D3 : L2M2NI 

However, not all the products of the model variables can be linearly independent: there 

are linear relations between: 

CiCj, AkDk, and AiAjEj, 

DiDj, AlB.Bj, and SfcC| (2.37) 

CiDi, AjBkCk, and A^BjCj 
for k a cyclic permutation of 1, 2, 3 and repeated indices are not summed. (It is plain 
that the three sets of products found to be linearly related must have the same Dynkin 
labels.) A specific form of the generating function, as expressed in terms of the elementary 
couplings, amounts to a specific choice of a set of forbidden couplings among those that 
are related by a linear relation. 
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3. Tensor-product descriptions 



3.1. The need for a tensor-product description 

It is clear that one major technical complication of the character method is that it 
starts at too fundamental a level, namely the character of the separate representations. 
One natural way to proceed is to start from a combinatorial description of the tensor- 
product rules. Such a description already takes into account the action of the Weyl group 
and encodes the various subtractions of the singular vectors. 

But how do we make the connection with the generating-function approach? The key 
is to find a combinatorial description which can be expressed as a set of linear Diophantine 
inequalities. Given this set of inequalities, there is an algorithm, again due to MacMahon, 
for constructing a generating function. (This is an adaptation of a method developed by 
Elliot [8] for the analysis of linear Diophantine equalities and for this reason the algorithm 
is often referred to as the EUiot-MacMahon method. For a detailed discussion of the 
algorithm, see in particular vol. 2 section VIII of [6].) This method is conceptually similar 
to the character method, except that the starting point is substantially closer to the end 
result. See section 7.3 for a slight generalisation of this algorithm. 

Although the description of tensor products via linear Diophantine equations is a more 
efficient route to finding the generating function than the character one, complications 
associated to the Q, projections remain a source of technical difficulty that severely limits 
the practical applicability of the method. 

A more powerful approach to our problem is to use the techniques of computational 
algebra. We start with a description of the tensor-product multiplicities as solutions to 
linear Diophantine inequalities. Efficient algorithms exists for finding the fundamental 
solutions to these inequalities [9] . From these we find directly a model for the generating 
function using Grobner basis techniques. (This is roughly the inverse of MacMahon's 
method which was originally conceived as a technique to generate the elementary couplings 
and their linear relations through the construction of the generating function. Here, the 
elementary couplings and their relations are first obtained and used as the input for the 
construction of the generating function.) 
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4. The LR rule (su{N)) 



For su{N) tensor products there is a particularly convenient description based on 
Littlewood-Richardson tableaux supplemented by the stretched-product operation (defined 
below) [10]. 

Integrable weights in su{N) can be represented by tableaux: the weight (Ai, A2, ■ ■ ■ , Ajv-i) 
is associated to a left justified tableau of A'" — 1 rows with Ai + A2 + ■ ■ ■ + Aat-i boxes in the 
first row, A2 + A3 + • • • + Aa^-i boxes in the second row, etc. Equivalently, the tableau has 
Ai columns of 1 box, A2 columns of 2 boxes, etc. The scalar representation has no boxes, 
or equivalently, any number of columns of N boxes. 

The Littlewood-Richardson rule is a simple combinatorial description of the tensor 
product of two su{N) representations A(S>/U. The second tableau (p) is filled with numbers 
as follows: the first row with I's, the second row with 2's, etc. All the boxes with a 1 are 
then added to the first tableau according to following restrictions: 

1) the resulting tableau must be regular: the number of boxes in a given row must be 
smaller or equal to the number of boxes in the row immediately above; 

2) the resulting tableau must not contain two boxes marked by 1 in the same column. 

All the boxes marked by a 2 are the added to the resulting tableaux according to the above 
two rules (with 1 is replaced by 2) and the further restriction: 

3) in counting from right to left and top to bottom, the number of I's must always be 
greater or equal to the number of 2's. 

The process is repeated with the boxes marked by a 3, 4, • • • , A?" — 1, with the additional 
rule that the number of i's must always be greater or equal to the number of i + I's 
when counted from right to left and top to bottom. The resulting Littlewood-Richardson 
(LR) tableaux are the Young tableaux of the irreducible representations occurring in the 
decomposition. 

These rules can be rephrased in an algebraic way as follows [10]. Define to be the 
number of boxes i that appear in the LR tableau in the row j. The LR conditions read 



k-i 



k 





(4.1) 
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and 

fc k 

^Ui-ij-i-^nij >0 2<i<k<N and i < N - 1. (4.2) 

The weight fi of the second tableau and the weight v of the resulting LR tableau are 
respectively given by 

AT N-1 

^'^ij=^l^j i = 1,2, ...,iV- 1 , 

(4.3) 

N-l min(j,iV-l) 
- ^j+^'^ij+l= ^ "^ij j = 1,2, ...,iV- 1 . 

Hence, given three weights A, ^ and z/, the number of non-negative integers solutions {riij} 
satisfying the above conditions gives the multiplicity M^^^ of in the tensor product A®//. 

The combined equations (4.1) and (4.2) constitute a set of linear and homogeneous 
inequalities. As described in [11], the Hilbert basis theorem guarantees that every solution 
can be expanded in terms of the elementary solutions of these inequalities. 

We can construct a model for the solutions of the equations (4.1) and (4.2) by intro- 
ducing new formal variables Ai, 1 < i < t where t is the total number of variables in (4.1) 
and (4.2). Then the subring of Q[Aj; 1 < i < t] generated by the monomials A" with a a 
solution of (4.1) and (4.2) provides the required model. This ring R will be generated by a 
finite set of monomials Ej 1 < j < s which we call elementary couplings corresponding to 
the elementary solutions of (4.1) and (4.2). Thus R is isomorphic to Q[ei, . . . , eg]// under 
the mapping (j) : ei ^ Ei where / is some ideal. Each element of / corresponds, via the 
map (/), to a relation between the elementary couplings. 

In the case of LR tableaux, there is a nice pictorial representation of the model R. 
Consider the set of formal linear combinations of LR tableaux with rational coefficients. 
It is given a ring structure by defining the stretched product of two LR tableaux (denoted 
by •) to be the tableau obtained by fusing the two tableaux and reordering the numbers 
in each row in increasing order [10]. More algebraically, if we denote the empty boxes of a 
LR tableau by a 0, so that 

N-l 

noj^J2^i j = l,2,...,Ar-l (4.4) 
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we can characterise completely a tableau by the data {n^} with now z > 0. It is clear 
the set of numbers {n^} with z > 0, or equivalently, {Aj,ny} with z > 1, is a complete 
set of variables for the description of the tensor products. Then, the tableau obtained by 
the stretched product of the tableaux {riij} and {n[j} is simply described by the numbers 
{uij + n'^j}- Here is a simple example: 







1 


1 


1 


2 


2 


3 




4 











1 




1 


2 




2 













1 


1 




1 


1 


1 


2 


2 




2 


2 


3 






A 













(4.5) 



This ring of tableaux is isomorphic to the model R constructed above and we do not 
distinguish between them. Thus we specify a set of elementary couplings (i.e. a set of 
generators of R) as a set of elementary LR Tableaux. 



4.1. Example: the su{2) case 

The complete set of inequalities for su{2) variables {Ai,nii, 77.12} is simply 

Ai > ni2 nil > ni2 > (4.6) 

The other weights are fixed by 

fJ'i = nil + ni2 1^1 = Xi+ nil - ^2 (4.7) 

By inspection, the elementary solutions of this set of inequalities are 

(Ai,nn,ni2) = (1,0,1), (1,0,0), (0,1,0) (4.8) 

which correspond respectively to £'i,£?2,-E'3 in (2.22). These correspond to the following 
LR tableaux: 



El 



E2: B, H 



(4.9) 



It is also manifest that there are no linear relations between these couplings. The generating 
function is thus simply: 



Qsu{2) _ 



il-Ei)(l-E2)il-Es) 
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(4.10) 



4-2. Example: multiple tensor products in the su{2) case 



Consider the problem of finding the multiphcity of the representation ^ in the triple 
product X <Si IJ, <Si I' D C- As a first step, the LR rule applies as before: with nn + 77-12 = 
III, we have Ai > 77-12. After the first product, we re-apply the LR rule with now Ai 
replaced by Ai + rin — n\2 and n^- replaced by 777.^^ with mn + 77112 = ^i- The LR gives 
Ai + 7111 — > 77712. The two inequalities for the su{2) quadruple product are then: 

Al > 72-12 Ai + ?7ii — 77i2 > 7T7i2 Uij > 777^^ > (4-11) 

The elementary solutions are then, in the order: name of the coupling, corresponding 
Dynkin labels and the 5-vector (Ai, 77ii, 7712, mn, 77712),: 

£^1 : (1) (1) (0) D (0) (1,0,1,0,0) 

£^2 : (1) ® (0) ® (1) D (0) (1,0,0,0,1) 

Es: (1)®(0)®(0)D(1) (1,0,0,0,0) 

(4.12) 

£^4 : (0) ® (1) (1) D (0) (0,1,0,0,1) 

£^5 : (0) ® (1) (0) D (1) (0,1,0,0,0) 

£^6 : (0) ® (0) ® (1) D (1) (0,0,0,1,0) 

The linear relation, whose existence was signalled by the character method, is 

£;3£^4 = £;2£^5 : (1,1,0,0,1), ^£;i£;6 : (1,0,1,1,0) (4.13) 

Choosing to forbid the product E2E4, the generating function can be written in the form 

1 — E3E4 



G = 



(1 - Ei)(l - E2){1 - Es)il - E4){1 - E5)il - Ee) 

(4.14) 



I n 1- Eil [1- Es^ 1- E4) 

yi=l,2,5,6 / 



The latter form makes manifest the absence of E3E4. 

We could represent the elementary couplings in terms of tableaux, where the boxes 
with I's refers to the tableau and those with 2's originate from the v tableau. (Warn- 
ing: the resulting tableaux describing the four- products are not necessarily LR tableaux.) 
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Hence, riij gives the number of I's in row j of the composed tableau while mik gives the 
number of 2's in row k. The elementary tableaux are 



El : 
E4: 



Eo 



,Es: □ 



1 
2 



From this representation, the relation reads 

-E3-E4 = E2E5 : 



Example: the stt(4) case 

The stt(4) LR conditions are: 

Ai > ni2 
A2 > ni3 

A2 + ni2 > ni3 + 77,23 

A3 > ni4 

A3 + nis > ni4 + n24 

A3 + ni3 + n23 > ni4. + ?^24 + ?^34 



7^ -Ei£^6 '• 



nil > n22 

nil + ni2 > 77,22 + n23 

nil + ni2 + 7113 > n22 + ?^23 + n24: 

n22 > 

n22 + n23 > n33 + 7134 



The tensor-product elementary couplings are: 

A^:H. ^3:0, Bi: 



Ai 



Ci:E, Q 



2 ■ 



C: 



3 : 



D' 



3 • 



together with 



: 





1 




1 






1 












2 










1 


2 


,1^2: 






,1^3: 








3^ 




1 






2^ 




2^ 





£^2 : 



E3 





1 




2 


1 




1 





(4.15) 



(4.16) 



(4.17) 



(4.18) 



(4.19) 
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The Dynkin-label transcription of the elementary couphngs reads 



Ai : 


(0,0,0) 


® (0,0,1) 


D 


(0,0,1) 


D[: 


(0,1,0) 


® (1,0,0) 


D (0,0,1; 


A2 : 


(0,0,1) 


(1,0,0) 


D 


(0,0,0) 




(1,0,0) 


® (1,0,0) 


D (0,1,0; 


A3 : 


(1,0,0) 


(8) (0,0,0) 


D 


(1,0,0) 


D'3: 


(1,0,0) 


® (0,1,0) 


D (0,0,1; 


Bi : 


(0,0,0) 


® (0,1,0) 


D 


(0,1,0) 


Di : 


(0,1,0) 


® (0,0,1) 


D (1,0,0; 


B2 : 


(0,1,0) 


® (0,1,0) 


D 


(0,0,0) 


D2 : 


(0,0,1) 


® (0,0,1) 


D (0,1,0; 


B3: 


(0,1,0) 


(0,0,0) 


D 


(0,1,0) 


D3: 


(0,0,1) 


(8) (0,1,0) 


D (1,0,0; 


Ci: 


(0,0,0) 


® (1,0,0) 


D 


(1,0,0) 


El : 


(1,0,1) 


8) (0,1,0) 


D (0,1,0) 


C2: 


(1,0,0) 


0(0,0,1) 




(0,0,0) 


E2 : 


(0,1,0) 


» (0,1,0) 


D (1,0,1) 


C3: 


(0,0,1) 


® (0,0,0) 


D 


(0,0,1) 


E3 : 


(0,1,0) 


» (1,0,1) 


D (0,1,0) 


(4) there are 


15 relations 


[12,10] : 














fc — CiEi 








EiEj — 


BkDkDf, 



(4.20) 



E>iEi = CjBf-Dk 
with i, j, k a cyclic permutation of 1, 2, 3 



D,Ei = BjD.Ck 



(4.21) 



To construct the generating function, we need to select forbidden couplings. It turns 
out that when there are more that one relation, complications may arise. We must ensure 
that the selected forbidden couplings are complete, which means that no further (usually 
higher-order) relations are required for a unique decomposition of a given coupling. A 
technique that is tailor-made for dealing with problems of that type is that of Grobner 
bases. This will be introduced in the next section. At this point, we simply indicate 
a complete choice of forbidden couplings, namely {E^Ej, D'^Ei, DiEi, DjD^, D'^Di}. This 
yields then a model for the generating function, which then reads [12,10] : 



i=l 



where 



+ D2^2^3^1 + Dibib2b3 + E3E3bib2 + D[b[biE3 

+ D'^E3b'^E3b'^ + E2E2b[b'^ + E2DiE2bib[ + E2D3E2b3b'^ 

+ DiDsE2biD3E2 + D2D'^D2b'^Ei + L'2^2^3^'2^'2^3)- 

Mi = {l-Mi)-\ 



(4.22) 



(4.23) 
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5. Diophantine inequalities: elementary couplings, relations and Grobner bases 



We introduce the idea of the Grobner basis via a simple example (see also [13]). Sup- 
pose i? is a model for a generating function, where R = Q[x, y, z, t]// and / = (xy—t, zy—t) 
is the ideal generated hy xy — t and zy — t, with an grading given by (1, 0), (0, 1), (1, 0) 
and (1, 1) for x, y, z and t. Writing x = x + I and similarly for the other variables, we have 
in R that xy = t and zy = t. These two expressions give two re-write rules : xy ^ t and 
zy H-^ t. These rules can be used to simplify any monomial. The aim is to find a re-write 
rule which, when iterated, produces unique representatives for the classes of /. If this is 
the case, then a vector space basis of R would consist of terms of the form m + I with m 
a monomial which is not divisible by any of the left-hand sides of the rewrite rules. 

In the example above, if we had 'good' rewrite rules then a basis for R would be 
represented by monomials not containing xy or zy, i.e. monomials of the form either y^'t^ 
or x°'z^t'^. The generating function which counts these monomials is: 

{\-AB) (t^^ (1-^)0 ' ^^'^^ 

The exponent of A carries the first grading index and B the second. 

However this generating function is not correct. It contains the term 2A?B corre- 
sponding to the monomials xt and zt. But the polynomial z{xy — t) — x{zy — t) = xt — zt 
is also in / and hence in R we have xi = zt and so the space of grade (2, 1) has dimension 
1 rather than 2. This problem can also be seen as a problem with the re-write rules. If we 
start with xyz then we can use the first re-write rule: xyz i— > tz or the second: xyz i— > xt. 
We cannot apply any further re-write rules and so this set of re-write rules does not pro- 
duce a unique representative. The solution is to include the rule xt zt. This gives a set 
of 3 rules: xy ^ t, zy ^ t and xt ^ zt. It turns out that this is a 'good' set and so a basis 
for R is given by (the classes of) monomials of the form y^'t^, x'^z^ and z'^t}' which gives 
the generating function: 

+ 7^+ . 1 ... (8.2) 



{\-AB){\-B) (1-^)2 (^\-A){\-AB) 

The set of 'good' generators, xy — t,zy — t, xt — zt we have found for / is known as a 
Grobner basis [14]. 
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The general procedure for constructing a Grobner basis given a set of generating 
polynomials is as follows. First choose a term ordering, which is an ordering on monomials 
with the property that any chain mi > m2 > . . . has finite length. For example we can 
order the variables hy x>y>z>t and then order all monomials by the corresponding 
lexicographic (dictionary) order, for example: x'^y > xyz > y^. For each generator of our 
ideal /, select the monomial which is highest with respect to the given term ordering. This 
is then the term which appears on the left of the re-write rule. The lexicographic ordering 
gives the first two re-write rules of our example: xy ^ t and zy t— > t. Next, for each pair 
of leading terms find the lowest common multiple and simplify it in the two possible ways. 
In this case there is only one pair of leading terms and the lowest common multiple is xyz 
which simplifies to xt and yt. Continue to apply the re- write rules until the terms do not 
simplify further. If the resulting pair of terms are the same, then proceed to the next pair 
of leading terms, otherwise add a new re- write rule. In this case we add xt ^ yt. Proceed 
until no pair of leading terms gives a new rule. This is the case for the rules we now have. 
For example the two rules xy ^ t and xt ^ zt appears to give a new rule by simplifying 
xyt to both and yzt. However the second term can be further reduced to and so no 
new rule is required. 

This algorithm for computing Grobner bases is known as Buchberger's [14] algorithm. 
Improvements on this basic algorithm mean that it is now feasible to find Grobner bases 
for quite large sets of generating polynomials. (The web pages of the computer-algebra 
information network at the address http: / /cand. can.nl/CAlN contain information about 



many of the programs currently available.) 

Although it is not clear from this example, Grobner bases are a very versatile tool for 
performing explicit calculations. We end this section with an illustrative example relevant 
to our discussion of tensor-product generating functions. 

Consider a set of linear Diophantine equations: 

Met = 0, a > (5.3) 

with M an integer matrix and a a vector of non-negative integers. We would like to 
construct a generating function for the solutions to this set of equations: 



(5.4) 



A non-trivial example is given by the Diophantine equations that describe a 3 x 3 magic 
square: 

(a b c\ 

(5.5) 

with non-negative entries and equal row and column sums. The magic square condition 
(the sum of each row and each column is the same, say equal to t) gives the following set 
of equations: 




a + b + c = t 
d+e+f=t 
g + h + i = t 



a + d + g = t 
b+e+h=t 
c+f+i=t 



(5.6) 



With a standing for the column vector with entries (a, b, c, d, e, /, g, h, i, t), the matrix M 
reads 



M = 



/I 


1 


1 




















-1\ 











1 


1 


1 











-1 




















1 


1 


1 


-1 


1 








1 








1 








-1 





1 








1 








1 





-1 


Vo 





1 








1 








1 


-1/ 



(5.7) 



There is a straightforward algorithm for finding the basic set of solutions [9] which 
yields: 

ai = (0, 0, 1, 0, 1, 0, 1, 0, 0, 1) a4 = (1, 0, 0, 0, 0, 1, 0, 1, 0, 1) 

a2 = (0,1,0,0,0,1,1,0,0,1) as = (0,1,0,1,0,0,0,0,1,1) (5.8) 

as = (0, 0, 1, 1, 0, 0, 0, 1, 0, 1) ae = (1, 0, 0, 0, 1, 0, 0, 0, 1, 1) 

We shall use A, B, . . . ,T to denote the "grading variables" of this example so that the 
exponent of A carries the value of a and so on. A model for the generating function is given 
by the subring S of Q[^, B, C, D, E, F, G, H, /, T] generated by monomials corresponding 
to the 6 elementary solutions. 

El = CEGT, E2 = BFGT, Es = CDHT, 

E4 = AFHT, E5 = BDIT, Eq = AEIT 

The monomials in S correspond to magic squares. For example EIE4^Eq = A^C'^E^FG'^HIT'^ G 
S corresponds to a square with row and column sums equal to 4: 



(5.10) 



(5.9) 




22 



Note that in this example it is convenient to construct our model as a subring of the 
ring of grading variables. Thus each "elementary coupling" Ei is actually equal to the 
corresponding monomial in the grading variables. 

However, there are relations between these generators and so it is not immediately 
clear how to construct the Poincare series for S. What we require is an isomorphism of S 
with R = Q[ei, . . . , eol/I such that Ei, i = 1, . . . ,Q and such that we have a Grobner 
basis of the ideal / (the 'ideal of relations'). 

Fortunately, such an isomorphism is easily constructed using Grobner-basis methods. 
Introduce the ring Q[^, B, C, D, E, F, G, H, /, T, ei, . . . , eg] with the lexicographic ordering 



Let J be the ideal generated hj Ei — ei, . . . , Eq — cq. This is not necessarily a Grobner 
basis with respect to this term ordering. Let G be the Grobner basis for J with the given 
ordering. Then it can be shown [14] that G fl Q[ei, . . . , ee] is a Grobner basis for the ideal 
of relations / which we require. In this case G is quite large, but its intersection with 
Q[ei, . . . , eg] is 616465 — 626366. The corresponding relation in R is EiE^Er, — E2E3EQ 
and these two terms do indeed give the same magic square, so that indeed we have found 
a relation between the generators of R. The Poincare series for Q[ei, ...,66]// is easily 
computed: 



6. Berenstein-Zelevinsky Triangles 
6.1. Generalities 

The previous examples make clear the usefulness of a re-expression of the tensor- 
product calculation in terms of Diophantine inequalities. The Littlewood-Richardson al- 
gorithm yields a set of such inequalities only for su{N). Fortunately, Berenstein and 
Zelevinsky [15] have expressed the solution of the multiplicity of a given tensor product 
as a counting problem for the number of integral points in a convex polytope. For a given 
algebra, the polytope is formulated in terms of a characteristic set of inequalities. For 
su{N), these reduce to the LR set of inequalities. For the other classical algebras, except 
sp(4), the proposed set of inequalities is a conjecture. 



A>B>C>D>E>F>G>H>I>T>ei>...>ee 



(5.11) 




(5.12) 
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6.2. BZ triangles for sp(4) 



The combinatorial description of tensor products for sp{4) is not as simple as in the 
su{N) case: a standard LR product must be supplemented by a division operation and 
modification rules [16]. Given the BZ set of inequalities, the natural way to proceed, as 
just mentioned, is to interpret these as the appropriate inequalities for the description of 
the tensor products. These inequalities are as follows: 



Ai >p 

A2 > ri/2 

A2 > ri/2 + q - p 

A2 > ^2/2 + q- p 

ui = r2 — ri — 2p + Xi + 



1^1 >Q 

Hi > q + ri - r2 
III >p + ri-r2 
1^2 > r2/2 

1^2 = p — q — r2 + X2 + fJ'2 



(6.1) 



(Our notation is different from that used in [15]: the relation is ri = mi, r2 — TO2, p = 
^12 5 Q = ''^'^12 •) The sp{4) tensor product coefficient Af\^i, is thus given by the number of 
solutions of the above system with ri, r2 G 2 N et p, g G N (N being the set of nonnegative 
integers). 

A proper set of variables for a complete description of a particular tensor-product 
coupling is thus {Ai, A2, A*2) i^i, 1^2, P, q}- We give the list of elementary couplings, adding 
to each coupling the corresponding four- vector [ri,r2,p,q]: 



Ai: 


(0,0) 


® (1,0) 


D 


(1,0) 


[0,0,0,0] 


A2: 


(1,0) 


® (0,0) 


D 


(1,0) 


[0,0,0,0] 


A3: 


(1,0) 


® (1,0) 


D 


(0,0) 


[0,0,1,1] 


Bi : 


(0,0) 


(0,1) 


D 


(0,1) 


[0,0,0,0] 


B2 : 


(0,1) 


® (0,0) 


D 


(0,1) 


[0,0,0,0] 


B3: 


(0,1) 


® (0,1) 




(0,0) 


[2,2,0,0] 


Ci : 


(0,1) 


(8) (1,0) 


D 


(1,0) 


[0,0,0,1] 


C2: 


(1,0) 


(8) (0, 1) 


D 


(1,0) 


[0,2,1,0] 


C3: 


(1,0) 


® (1,0) 


D 


(0,1) 


[0,0,1,0] 


Di 


(2,0) 


® (0,1) 


D 


(0,1) 


[0,2,2,0] 


D2 


(0,1) 


«) (2,0) 


D 


(0,1) 


[2,0,0,0] 


D3 


(0,1) 


(0, 1) 


D 


(2,0) 


[0,2,0,0] 
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(6.2) 



The unspecified finear relations mentioned in (2.37) can now be obtained. To find those 
products that are equal in the current situation we need only compare their corresponding 
sets of four- vectors [ri, r2,p, q\ (which are additive in products of couplings). We thus find 
for instance that 

C1C2 = A^D^ : [0, 2, 1, 1] ^ ^1^253 : [2, 2, 0, 0] (6.3) 
Proceeding in this way for the other cases, we find the following complete list of relations: 

C1C2 = A^D^, C2C3 = ^iDi C3C1 = ^1^352 

DiD2 = BsCi D2Ds = AlB2Bs DiD^ = B^Cl (6.4) 
Ci-Di = A3S2C2 C2-D2 = AiB^C^ C^Ds = A1B2C2 

The use of the BZ inequalities to find the elementary couplings and their relations is 
novel. (An off-shoot of our construction is that it provides an indirect proof of the validity 
of the BZ inequalities since we recover from it the result of [7] derived from the character 
method.) 

A possible choice of forbidden products is the one given in [7]: 

{CiC^,DiD^,CiDi} (6.5) 
with i, J = 1, 2, 3 and i ^ j- It leads to the generating function: 

\ti J ^ (6.6) 



+C2C2D3 + D1C3D1 + C3C3D2) 



Of course, by modifying the ordering in the Grobner basis, we can get other choices of 
forbidden couplings. Here is another set of forbidden couplings that can be obtained: 
{DiDj,CiDi, AiDi, A3D3, AiA^B2}. The corresponding generating function reads 



QSp{4) ^ B1B2B3 



W aA (7,(1 - A1A3S2) + D3D3A1A2C1C2 



(6.7) 

^Dxb^A2A3C2C3 + D2b2hA2MCxC3{\ - A1A3S2) . 

These two generating functions are equivalent when rewritten in terms of the grading 
variables, that is, in terms of Dynkin labels. However, they originate from two distinct 
models. The second one turns out to be well adapted to the fusion extension. 
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7. A vector basis approach to the construction of generating functions 

In this section, we present a simple and systematic way of generating by hand all the 
elementary solutions of a set of linear homogeneous inequalities starting from the well- 
known construction of a vector basis. The first step amounts to reformulate the system of 
inequalities in terms of equalities. We then look for the elementary independent solutions 
by relaxing the positivity requirement. In other words, we construct the vector basis. In 
a final step, we find the minimal linear combinations of these vector basis elements that 
yield positive solutions. This will also provide an illustration of MacMahon's projection 
technique. The result of this projection is the desired tensor-product generating function. 
Hence, this approach turns out to be a new way of constructing the tensor-product gen- 
erating functions. (This generic method, referred to as being novel for tensor products, is 
certainly well-known in general: it is discussed in the first reference of [11].) 



7.1. Graphical representations as BZ triangles for su{N) 

Consider the direct transformation of the LR inequalities to equalities by introducing 
an appropriate number of new non-negative integer variables. Consider first the su{2) 
case, for which there is a single inequality: Ai > ni2. We transform this into an equality 
by introducing the non-negative integer a defined by 

Ai = ni2 -I- a (7.1) 

The expression for v\ becomes then i^i = Ai -|- nn — n\2 = a + nn. Since ni = un + n\2i 
we are led naturally to a triangle representation of the tensor product: 

X<^uDv ^ (7.2) 
ni2 nil 

We read ofi^ the Dynkin label of the A representation from the sum of the two integers 
that form the left side of the triangle, that of the // representation from the bottom of the 
triangle and the ui label is the sum of the two integers that form the right side. A more 
uniform notation amounts to setting a — mi2 and nn = /12, in terms of which the triangle 
looks quite symmetrical: 

(7.3) 

ni2 I12 
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with 

Ai = mi2 + ni2 111 = ni2 + I12 J^i = toi2 + h2 (7.4) 

These numbers mi2 and /12 plays the role of ni2 in the permuted versions of the tensor 
product. The triangle combinatorial reformulation of the tensor product problem is as 
follows: the number of triangles that can be formed from non- negative integers ni2, mi2 
and li2 that add up to the Dynkin labels of the representations under study according to 
the above relations gives the multiplicity of the triple coupling A D i^, or equivalently, 
the multiplicity of the scalar representation in the product A®//®^' D (0) (since for stt(2), 
ly* = v). 

The situation for sm(3) is somewhat more complicated. The transformation of the LR 
inequalities (4.1, 4.2) into equalities in this case takes the form 

Ai = ni2 + a nil = n22 + d 

A2 = ni3 + 6 nil + ^12 = ^22 + ri22, + e (7.5) 
A2 + ?^l2 = rii^ + ?i23 + c 
The expression for the other weights becomes 

A*! = ni3 + e 11,2 = n22 + "-23 

(7.6) 

1^1 = a + a V2 = n22 + c 

Since there are two expressions for both nn and A2, there follows the compatibility rela- 
tions: 

^^12 -\-d = n23 + e n23 + c = b + ni2 (7.7) 
By adding these two relations, we find: 

c + d=b + e (7.8) 

Again we are led naturally to a triangle representation: with C = ^* this reads 



a 

ni2 d 

b c 
ni3 e n23 ^22 



(7.9) 



We read the Dynkin labels from the sides of the triangles, from Ai to C2 in an anti-clockwise 
rotation starting from the top of the triangle, exactly as for su{2), except that here there 
are two labels on each sides. Notice that the compatibility conditions amounts to the 
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equality of the sums of the extremal points of the three pairs of opposite sides of the 
hexagon obtained by dropping the three corners of the triangle. 

Again a more symmetrical notation is: 

a = mi3 6 = 777,23 C = mi2 d = l23 e = li2 ^22 = ^13 (7-10) 

in terms of which the triangle reads 



mi3 

"7,23 1^12 
"13 h2 "23 ^13 



(7.11) 



with labels fixed by: 



Al = mi3 + 7712 A2 = 777,23 + "l3 
Atl= "13 + ^12 A*2 = "23 + ^13 (7-12) 
Cl = ^13 + 1^112 C2 = ^23 + "il3 

The hexagon conditions read: 

"12 + TO23 = "23 + mi2, 

^12 + "i23 = ^23 + "^12, (7-13) 
^12 + "23 = ^23 + "12- 

In terms of triangles, the problem of finding the multiplicity of the su{3) tensor product 
A (8) (8) C ^ boils down to enumerating the number of triangles made with non-negative 
integers that form a bipartition of the Dynkin labels and that satisfy the above three 
hexagon relations. (For the su{N) generalisation, see [17]). 

Here is the rationale for the labelling nij^rriij, lij from the triangle point of view [18]. 
If Ci are orthonormal vectors in R-^, then the positive roots of su{N) can be represented 
in the form Cj — e-,-, I < i < j < N. The triangle encodes three sums of positive roots: 

li + C-X* = ^ kj{et-ej) , 

i<j 

C + X-fi* = mij{ei-ej) , ^^^^4^ 

i<j 

X + li-C* = ^ nij{ei-ej) , 

i<j 

The hexagon relations are simply the consistency conditions for these three expansions. 
Clearly, the variables Uij that appear in the above relations are exactly the Uij that appear 
in the LR tableaux for the product A (8) D C* = 
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1.2. From a vector basis to the generating function: the su{3) case 



Given the transcription of inequalities into equalities, we can easily extract the cor- 
responding basis vectors. This is the starting point of a new method for constructing the 
tensor-product generating functions. To keep things concrete, we focus on the su{3) case. 
The goal is to first get a vector basis and then to project it to get the elementary couplings. 
The generating function is a direct result of this procedure. 

The equality version of the LR inequalities are (7.12) and (7.13); they underlie the 
construction of the BZ triangle (7.11). The last hexagon condition of (7.13) is the difference 
of the previous two so it is not an independent relations. We thus have a total of 15 
variables: Ai, • • • , ^2, li2, • ■ ■ , ^^23 and 8 equations. The number of independent variables is 
thus 7. These will be chosen to be mi3, m23, hs, I23, ^^12, nis, 77-23. The dependent variables 
are fixed as follows: 

Ai = mi3 -I- ni2 A2 = m23 + ^13 

A*! = nis + ni2 + I23 - n23 A*2 = ^23 + ^i3 

Ci = "-12 + + ha - n2s C2 = ^23 + "7,13 

h2 = ni2 + I23 - n23 rni2 = ni2 + m23 - 'n'23 



(7.15) 



We now look for the elementary solutions of this system (without invoking the constraint 
that all the above dependent variables should be necessarily positive). The sought basis 
vectors are obtained by setting one of the variable mi3, • • • , 77-23 to 1 and all other set equal 
to zero. 

This produces (in order) the triangles £"2, £"5, ^e, -E's, -E"?, -E'4 and Zi displayed below: 

£2:(1,0)(0,0)(0,1) £3:(0,0)(1,0)(0,1) 

1 
1 



1 

£4 : (0,1)(1,0)(0,0) £5 : (0, 1)(0, 0)(1, 0) : (0, 0)(0, 1)(1, 0) 





11 

1000 0000 0001 
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E, : (1,0)(1,0)(1,0) Zi : (0, 0)(-l, 1)(-1, 0) 

^ ^ 

1 ('^•16) 

1 0-1 
1 - 1 1 

These are all genuine BZ triangles except for Zi which has some negative entries. However, 

at this level, there are no relations between these elementary solutions (the basis vectors 

are independent), hence the decomposition of any solution in terms of these 7 basic ones 

is unique. All solutions are then freely generated from the following function: 

^ = (1 - E2){1 - E^){1 - E^){\ - E^){\ - E^){\ - E^){\ - Zi) ^^-^^^ 

To recover the generating function for all tensor products from the above expression, we 
need to project out terms that lead to triangles with negative entries. To achieve this, 
we introduce the grading variables associated to the above couplings (compare the above 
triangles with the general form given in (7.11)): 

E2 '■ Mi3, £^3 : L12L23 E4 : A^i3 

E5 : M12M23 Ee : L13 E^ : L12M12N12 (7.18) 

Zi : L-^^M-^N23 

Our generating function follows from the projection of the above function G, re-expressed 
in terms of the grading variables, to positive powers of L12 and M12. Equivalently, one 
can re-scale L12 by x and M12 by y and project to positive powers of x and y and set 
X — y — 1 in the result. This is equivalent to the rescaling 

Es^xEs E5^yE5 Er ^ xyEr Zi^x'^y'^Z^ (7.19) 

We are thus led to consider 

h!lG{E2,xE3,---,x-^y-^Zi) (7.20) 



> > 



Keeping only those terms which depend explicitly upon a; or y, we have then 

X y 1 



>> (1 - xE3){l - yE^){l - xyEj){l - x'^y'^Z^) 



1 / 1 x-^y-^Zi 



+ 



{1 - xE3){l - yEr,){l - E^Zi) \l-xyE7 1-x-^y-^Zi 
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(7.21) 



No more work is needed for the first term. For the second one, we have 



X y 



X ^Zi f yE^ 1 



+ 



>> {l-xE3){l-ErZi){l-x-^ZiE5) \l-yE5 1-x-^y-^Zi 

x-^E^Zi 



X 

= Q 



> (1 - E,){1 - E7Zi){l - xEs){l - x-^ZiE^) 
x-^E^Zi ( xEz 



X 



+ 



> {l-E^){l-EjZ^){l-E3E5Zi) \l-xE3 l-x-^ZxE^ 
_ E3E5Z1 
~ (1 - £^5)(1 - ErZi)(l - EsE^ZiXl - E3) 

We then introduce the following two new elementary couplings 

El = EfZi , Eg — E^E^Zi 

Collecting the two terms resulting from the projection, we end up with 



(7.22) 



(7.23) 



(7.24) 



which is indeed the s«(3) tensor-product generating function. 



It is worth pointing out that the EUiot-MacMahon algorithm that has been presented 
here as a method distinct from the vector basis, can be reinterpreted in a way that makes 
the two approaches equivalent. This is done in section 3 of the first reference of [11]. There, 
the elementary solutions are not obtained as above by setting successively one dependent 
variables equal to 1 and the others equal to 0, but in reading them off directly from the 
columns of the 8x7 matrix of the matrix version of the above equation: 



/I 


1 






















1 


1 














1 





1 


-1 





1 














1 


1 








1 


1 





-1 


1 





1 

















1 





1 








-1 





1 


Vo 


1 


1 





-1 





0/ 



ni2 
ni3 

V hz ) 



( \x\ 

A2 

Ci 

C2 

mi2 



\mi2 / 



The exponentiated version of the columns gives the elementary solutions written below. 
This leads to the so-called 'crude' generating function that is then projected onto the 
positive solutions by the usual method. 
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1.3. General aspects of the vector basis construction 



In general, of course, the fundamental solutions to the linear system may have non- 
integral values of the variables. However the corresponding terms in the generating func- 
tion can be eliminated by rationalising all the denominator terms and then keeping only 
those terms in the numerator that have integral exponents. This suggests the following 
modification of MacMahon's algorithm. 

Consider the system of equations 

Ma; = 0, xeN'' (7.25) 

where M is a matrix of rank s. We thus have k variables and s relations between them. 
The dimension of the vector basis is thus k — s. We will denote the independent (free) 
variables as Xi, i = 1, ■ ■ ■ ,k — s and the remaining ones as Xj, j = 1, ■ ■ ■ , s. To find a 
generating function for the solutions of this system: 

1. First construct a basis in Q'^ for the solutions of Mx = by setting Xi = 1 with all 
other Xj zero (j = 1, ■ ■ ■ , — s, j 7^ i). Denote by x^^"* the value of the dependent 
variable Xj evaluated aX Xi — 1 with all other Xi zero. The basis then reads 

ei = (l,0,0...,0;{£5')}), 
e2 = (0,l,0...,0;{xf }), 

(7.26) 

efe_, = (0,0,0,...,l;{£f-^)}) 

By construction, the e^'s are linearly independent. However notice that in general the 
x^j^ might be rational. 

2. From the form of the e^'s, it follows that any solution to (7.25) can be written as 

CiSi with Ci non- negative integers. In particular this means that every solution to 
(7.25) corresponds to a term in the generating function: 

^^^^ = il-X^^){l-X^^)...{l-X^s) (^-27) 
where Xi, . . . , Xj^ are grading variables. 

3. G{x) may contain negative or fractional exponents due to the occurrence of x^^^ in the 
exponents. These are eliminated by first using MacMahon's algorithm to eliminate 
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any negative exponents and then rationalising denominators and keeping only terms 
with integral exponents in the numerators. 

The result is the generating function for the solutions to (7.25). This algorithm, 
however, does not seem to be optimal in all case. 

l.J^.. Multiple su{2) products from the vector basis construction 

A simple and different application of the formalism just developed is furnished by the 
analysis of su{2) quadruple tensor products. This application is different in that it does 
not rely on the triangle description and as such, its formulation is less direct. This does not 
mean however that there are no diagrammatic representations for the quadruple product. 
In fact, having a set of inequalities, wc can transform then into equalities, as it is done 
below, and from them set up a diagrammatic representation. In the present case, it would 
correspond to two adjacent su{2) triangles, one upside down, with their adjacent sides 
forced to be equal. However, our analysis will not rely on such a description. It will serve 
as a preparation the somewhat more complicated sp(4) example treated in the following 
section. 

The Diophantine description of this problem has been presented in section 4.2. It is 
based on the two inequalities (4.11) which are readily transformed into equalities by the 
introduction of two non-negative integers ai, a^'- 

Ai = ni2 + ai Ai + rin - ni2 = mi2 + 02 (7.28) 

However this system does not contain any reference to the variable mn and for this reason 
we introduce the further constraint mn > which calls for a new non- negative integer 
variable: 

mil = 03 (7.29) 

We have thus a total of 8 variables : {Ai, nn, 77-12, mn, mi2, ai, a2, 03} and 3 equations. 
There are thus 5 independent variables, chosen to be {01,02,03,77,12,777-12}. The basis 
vectors, with components ordered as follows 

(ai, a2, as, ni2, 77712; Ai, rin, rrin) (7.30) 
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are obtained by successively setting equal to 1 one of {01,02,03,7112,777,12} and the others 
equal to 0. These basis vectors together with their exponentiated version written in terms 
of appropriate grading variables read: 

(1,0, 0,0,0; 1,-1,0) :LiAr-Ui 

(0,1, 0,0,0; 0,1,0) :iVii^2 

(0,0, 1,0,0; 0,0,1) :Mn^3 (7-31) 

(0,0, 0,1,0; 1,0,0) :LiiVi2 

(0,0, 0,0,1; 0,1,0) :iViiMi2 

The desired generating function is obtained from the projection to positive powers of A^n 
of the function 

1 

(J 32) 

(1 - LiiVfiUi)(l - iVii^2)(l - i^iiVi2)(l - iViiMi2)(l - Mii^s) 

The projection operation is done by the familiar method and the result, after setting all 
A = 1 is 

G= 1 - LiiVnMi2 . . 

(1 - LiiVi2)(l - LiMi2)(l - Li)(l - iVnMi2)(l - iVn)(l - Mn) ^ ' ^ 

from which we read of the 6 elementary couplings Ei, ■ ■ ■ , Eq (in the order where they 
appear in the denominator) given in (4.12) and the relation EsE4 = E2E5. The above 
function is exactly the one derived in section 4.2. 

7.5. sp{4:) diamonds and the vector basis derivation of the generating function 

The system of inequalities (6.1) pertaining to sp{4) can be transformed into a system 
of equations in the standard way: by setting ri/2 = Si and r2/2 = S2 and introducing the 
non-negative integers Oj, we get [19]: 

Ai = p + Oi ZV2 = 04 + Og 

A2 = si + 02 a2+p = a3+q 

//i = g + 05 03 + 81=04-1-82 (7.34) 

At2 = 52 + og 05 -I- 2S2 = O6 + 2si 

ui = ai+ar a6 + q = ar +p 
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This leads to a diamond-type graphical representation of the tensor product that has the 
advantage over the one presented in [15] of being linear in that the sum of two diamonds 
is also a diamond. This is illustrated in Fig. 1. 

In Fig. 1, all data pertaining to the first (second) Dynkin label appear at the left 
(right). Dotted lines relate those two points that compose the label indicated beside it. 
Opposite continuous lines are constrained to be equal, with the length of a line being defined 
as the sum of its extremal points except for the lines delimited by the points (ag, si) and 
(as, S2) where the point Si is counted twice (the little bar besides si and .s'2 being a reminder 
of this). Explicitly, for those lines, we have thus the constraint ae + 2si = 05 + 2s2- Given 
a triple sp{4) product, the number of such diamonds that can be drawn with non- negative 
entries yields the multiplicity of the product. For instance, the two diamonds that describe 
the triple coupling (1, 1) (1, 1) ® (2, 0) are shown in Fig 2. 

The dimension of the vector basis is 8 (18 variables and 10 equations, the last four 
equations above being linearly independent). As our free variables we choose the set 
{si, S2,p, q, ai, as, ae, ag}. The 8 basis vectors in terms of grading variables are: 

\2 o c . j\/r— Sji/T AT— 1 A—1 A— 2, 



£i : L2MtN2A4AiSi £2 : M^^M2N^'A2'A^^S2 

£3 : LiL^^A^f M^Uf £4 : L2M1N1A2A7Q 

£5 : LiNiAi £5 : L2N2A2A3A4 

£7 : MiNiA5A6A7 £8 ■■ M2N2A8 



(7.35) 



The generating function is obtained by first projecting of the function ]^(1 — £i)~^ to 
positive powers for each grading variables and then by setting all grading variables equal 
to 1 except for Li, Mi, NiS. The sp(4) elementary couplings are simple products of the 
£iS (the following ^1,2,3 should not be confused with the above grading variables): 



Ai 


= £7 


A2 


^£5 


As 


= £s£a 


Bi 


= £s 


B2 


= £6 


Bs 


= £i£2 


Ci 


= £4 


C2 


= £2^3^6^7 


C3 


= £i£z£7 


Di 


= £2£i£i£7 


D2 


= £1 




— £2£q£7 



(7.36) 



The complete list of sp{A) elementary couplings (6.2) are thus recovered. 
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8. Conclusion 



As was stressed in the introduction, the main purpose of this work is to prepare the 
ground for the analysis of fusion rules, which is the subject of a sequel paper. In this 
paper, we have reviewed the existing techniques for computing tensor-product generating 
functions and presented a comparative assessment of their virtues and limitations. We 
also focused on a model formulation linking generating functions to Poincare series, an 
idea first introduced in [11] and extended in [10]. Our contribution has been to rephrase 
this program more explicitly, clarify some issues and to exemplify the procedure with many 
examples, some of which are new. An extended version of this article is available on the 
Los Alamos server [20] . 

Acknowledgement: We thank R.T Sharp, J. Patera and M. Walton for useful discus- 
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